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1. INTRODUCTION 



In spite of the impressive progress that Cosmology has experienced in recent years, we 
are still missing a consistent explanation of the origin of the Universe and the formation 
of structures which is deduced entirely from a fundamental theory. General Relativity 
(GR) breaks down in the very initial instants of the history of the Universe, leading to a 
cosmological singularity of the big bang type I n this regime GR cannot be trusted, 
and the very own predictability of the laws of physics is lost. One expects instead that 
the physics of the Primitive Universe belongs to the realm of Quantum Gravity, namely, 
a theory of the gravitational field which incorporates the quantum behavior of nature. 
One of the most promising candidates for such a theory is Loop Quantum Gravity 
At present, important efforts are being made in order to adapt the techniques of Loop 
Quantum Gravity to much simpler settings than those of the complete theory, which on 
the other hand remains to be concluded. This is the case of a series of cosmological models 
obtained from GR by symmetry reduction. The resultingfield of research is known under 
the general name of Loop Quantum Cosmology (LQC) 3]. 

The first cosmological system whose quantization was performed to completion in 
LQC was the homogeneous, isotropic, and spatially flat model provided with a minimally 
coupled, homogeneous, and massless scalar field as matter content 
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. The geometry 

was polymerically quantized (i.e., using typical LQC methods), while the matter field was 
described using standard quantization methods. A thorough analysis of the resulting 
quantum dynamics 

aa 

showed that the initial singularity is successfully resolved. The 
classical big bang is replaced with a quantum big bounce which deterministically connects 
a semiclassical expanding universe with a previous semiclassical contracting one. After 
the pioneer study of this simple model, LQC has been further developed to describe other 
homogeneous isotropic systems (7f, or homogenous and anisotropic cosmologies like, e.g., 
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the Bianchi I model [8|, 121, HOj. In all these works, the cosmological singularity is eluded in 
the quantum dynamics, a result which strengthens the validity of the singularity resolution 
mechanism and the relevance of LQC. 

One may wonder whether the quantum resolution of the singularities of GR is just 
an artifact of the great symmetry of the models analyzed so far in the literature, and in 
particular of the homogeneity, or whether, on the contrary, singularities are removed as 
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well when inhomogeneities are present. In order to ask this question, it seems unavoidable 
to extend the quantization performed in LQC to inhomogeneous scenarios. Moreover, 
since inhomogeneities may have played a crucial role at the first instants of the Universe, 
one should take them into account, indeed, in the development of a realistic theory of 
quantum cosmology. 

In order to progress in this direction, we will consider here the quantization of one of 
the simplest inhomogeneous cosmological systems, namely the linearly polarized Gowdy 



T 3 model 
On the one 



111 ]. This model is a natural test bed to incorporate inhomogeneities in LQC. 



rand, its quantization by means of standard techniques has been discussed 



in detail 12], and a successful Fock (and Schrodinger) quantization has already been 



achieved 



141 ] . Even though this system has no timelike isometry and possesses an 



infinite number of degrees of freedom, this quantization has been shown to be unique (up 
to unitary equivalence) under certain reasonable requirements, namely, it is the unique 
Fock quantization with a unitary implementation of the dynamics and invariance under 
the action of diffeomorphisms (strictly speaking, of the only diffeomorphisms which remain 
in the system after a natural gauge fixing: the S l translations). On the other hand, 
generically, the classical solutions of the Gowdy model represent spacetimes with an initial 



curvature singularity 
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161 ] . Besides, the subset formed by the homogeneous solutions 



represents Bianchi I spacetimes with a three-torus topology, a model which has been 
polymerically quantized and where the singularity has again been shown to be resolved as 
a result of the quantization qI, lo|. Therefore, it is natural to ask how the inclusion of the 
inhomogeneities affects the singularity resolution in this particular model. In view of the 
knowledge available about the system and the characteristics of the different gravitational 
degrees of freedom entering the model, the simplest possibility to investigate that question 
is to perform a hybrid quantization, which combines the polymeric quantization of the 
homogeneous degrees of freedom with the Fock quantization of the inhomogeneities. This 
idea constitutes the basis of our work, whose main aspects have already been presented 
in Ref. 17J. Actually, this hybrid quantization will allow us to investigate not only 
whether the singularity is resolved owing to quantum geometry effects, but, furthermore, 
to explore whether the loop quantization of just the degrees of freedom that parametrize 
the homogeneous solutions (zero modes in our description) may suffice to cure the big 
bang singularity. As we will see, the answer turns out to be in the affirmative for our 
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model. The simplicity of this result reinforces the interest of our hybrid approach, since 
the same strategy can be applied to the quantization of more general cosmologies. 

In addition, the hybrid quantization of the Gowdy model provides a specially suitable 
arena for the analysis of other important issues in quantum gravity and cosmology, such 
as: 

• The role of the internal time. We will select as emergent time a gravitational 



theory. However, in contrast with 
we will quantize this emergent 



variable which behaves like a clock in the classical 
the situation studied in the recent literature 0, ^], 
time using a polymeric representation, because it describes one of the degrees of 
freedom of the homogeneous sector of solutions. As a consequence of this polymeric 
quantization, we will see that, in principle, the evolution with respect to that time 
is not unitary. 

• The recovery of the standard quantum field theory. One can regard the inhomo- 
geneities as (conveniently scaled) gravitational waves propagating in a homogeneous 
curved background, namely a Bianchi I spacetime. In a standard quantization, the 
inhomogeneities can be described by means of a Fock space. One may then ask 
what is the status of this Fock description from the viewpoint of LQC, where one 
adopts a polymeric quantization inequivalent to the standard one. We will see 
that, with our hybrid approach, one indeed recovers the standard Fock description 
of the inhomogeneities, even though the background is quantized employing loop 
techniques. 

The rest of the paper is organized as follows. In Sec. 2 we will introduce the classical 
system. We will partially fix the gauge and parametrize the degrees of freedom of the 
homogeneous solutions (what we will call the homogeneous sector) in Ashtekar variables, 
whereas the inhomogeneities will be described by nonzero Fourier modes of a certain scalar 
field (this will be called the inhomogenous sector). A global diffeomorphism constraint 
remains in the reduced model, affecting only the inhomogeneous sector. This sector will 
be quantized in Sec. 3 adopting a Fock representation. In addition, we will impose the 
quantum analog of the mentioned diffeomorphism constraint. On the other hand, a global 
Hamiltonian constraint, coupling the homogeneous and inhomogeneous sectors, is also 
present. In Sec. 4 we will construct a polymeric representation of the homogeneous sector, 



which coincides with the phase space of a vacuum Bianchi I model. The Hamiltonian 
constraint of the complete system will be imposed in Sec. 5, where we will also construct 
the physical Hilbert space of the hybrid Gowdy model. Finally, in Sec. 6 we will discuss 
the main results of our quantization. 



2. THE CLASSICAL GOWDY MODEL 



The Gowdy cosmologies are vacuum spacetimes with two spacelike commuting Killing 
vector fields and spatial sections of compact topology 111]. We consider the simplest of 
these Gowdy models, namely, the case with the spatial topology of a three-torus, T 3 , 
and with a linear polarization. In this case both Killing vectors are axial (because of the 
topology) and hypersurface orthogonal (owing to the considered polarization). 

Let us denote by d a and dg the two Killing vector fields. To describe the system, we 
choose global coordinates {t, 9, cr, 5} adapted to the symmetries, with 8, a, 5 G S 1 . The 
metric components depend on t and 9, being periodic in the latter coordinate. The space- 
time is globally hyperbolic and, in a 3+1 decomposition, the metric can be described by 
the densitized lapse function N, the shift vector N l , and the three-dimensional metric 
qij induced on the spatial slices which foliate the manifold, where i,j G {9, a, 5}. As a 
consequence of the hypersurface orthogonality, the induced metric is such that q a $ = 0. 
Besides, the conditions that q$ a and qes vanish fix the gauge freedom associated with the 
diffeomorphism constraints in a and S, the directions defined by the Killing vector fields 
[hs! . In addition, the dynamical stability of these conditions imply that N a and N 5 must 
be equal to zero. At this stage, the induced metric is diagonal and therefore can be char- 
acterized by three fields that describe the norm of one of the Killing vectors, the area of 
the isometry group orbits, and the scale factor of the metric induced on the set of group 
orbits. We further demand that both the generator of the conformal transformations of 
this latter metric and the area of the isometry group orbits be homogeneous functions. 
These conditions turn out to fix the gauge freedom associated with the nonzero (inhomo- 
geneous) Fourier modes (with respect to the ^-dependence) of both the diffeomorphism 
constraint in the ^-direction and the densitized Hamiltonian constraint. Besides, they 
imply that N e and N must be homogeneous functions Q]- The difference between 



this gauge fixing and reduction procedure and the one considered in Ref. is that 
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in that reference, the system was totally deparametrized, since the area of the isometry 
group orbits (which is a global time function) was entirely fixed, whereas now we leave 
unfixed its homogeneous part (the zero Fourier mode). 

As a result of the almost complete gauge fixing that we have performed, all the gauge 
degrees of freedom are fixed except two. The resulting reduced phase space splits into two 
sectors, the homogeneous sector and the inhomogeneous one. The former of these sectors 
is formed by the degrees of freedom that describe homogeneous metric functions in our 
gauge fixing, together with their conjugate variables. This includes the zero mode of the 
only metric field which has not been fixed yet, and which provides the norm of one of the 
Killing vector fields (e.g., dg). The inhomogeneous sector, on the other hand, contains the 
information about all the nonzero modes of this metric field and their conjugate momenta. 

We next describe the homogeneous sector in terms of Ashtekar variables: an SU{2) 
gravitational connection A\ and a densitized triad E l a , both of them constant on spatial 
sections. It is worth commenting that this homogeneous sector can be interpreted as the 
phase space of a vacuum Bianchi I model whose compact sections are homeomorphic to 
a three-torus. These are precisely the spacetimes described by the considered degrees of 
freedom when the inhomogeneities vanish. Hence, the system is provided with a natural 
coordinate cell, namely the T 3 -cell, with sides of coordinate length equal to lix . In a 



diagonal [SU (2)] gauge, and choosing from now on the Euclidean metric as fiducia 



the nontrivial components of the corresponding Ashtekar variables 9|, |2Cj are 21] 



metric, 



At = ^f, K = £,51 (1) 

The nonvanishing Poisson brackets between these variables are {c l ,pj} = 8irGj8j. Here, 
a is an internal SU (2) index, the symbol G stands for the Newton constant, and 7 is the 



Immirzi parameter 



22]. 



The rest of degrees of freedoms of our reduced system correspond to a metric field 
£(#) and its conjugate momentum P^(6) 13], both of them devoid of any zero mode con- 
tribution since these have already been included in the homogenous sector. Appropriate 
variables for the description of what we call the inhomogeneous sector are then all the 
nonzero Fourier modes {(^ m ,P™),m e Z — {0}}. Now, we introduce the creation and 
annihilation variables {(a m ,a* m )} which would be naturally associated with if this 
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were a free massless scalar field. These are defined by 



7T 



\G\m\ 



TV 



and the complex conjugate relation, and are such that {a m , a*^} = —iS mm- 
The explicit relation between the chosen variables and the metric is fla . 



ds 2 ^ \P0P«PS\ 



4tt 2 
where 



N 2 i2 d9 2 \ -^da 2 -%i-d5' i 
(27r) 4 / p 2 p 2 i 



and 



with 



= E -^ a ™ + O^, (4) 

m^O V I I 

= , _2<^ - f>)] 2 - -^^CW, (5) 

\c a p a + c s p 5 J J\pe\ \Pe\ c a p a + c s ps 



C( e ) = i E si g n ( m + m ) + m — ( a -m ~ <4) (a m +m + a- (m+ ^))e' me . (6) 



If the inhomogeneities vanish, we recover the Bianchi I spacetime. Let us comment that 
7 contains a zero mode contribution arising from quadratic terms in the inhomogeneities. 



We note also that pg is proportional to the time variable chosen in Ref. 13] to de- 
parametrize the system. 

The gauge fixing that we have performed is partial. This explains the appearance of a 
homogenous densitized lapse N in the metric, and the fact that is still allowed to redefine 9 
by introducing a homogeneous shift N e . As a consequence, two global constraints remain 
in the model. One of them is the zero Fourier mode of the diffeomorphism constraint 
in 9, Cq, which generates translations in the circle. The other is the zero mode of the 
densitized Hamiltonian constraint, Cq. The volume of the Bianchi I counterpart of the 



Gowdy spacetime (for vanishing inhomogeneities), which is given by V = \f\peP~aPs\ and 
that we will call the "homogeneous- volume" , allows us to define a global (nondensitized) 
Hamiltonian constraint, Cq, instead of Cq. It is obtained by choosing the (nondensitized) 
homogeneous lapse function N = VN rather than N. The two remaining constraints can 
then be expressed as \v\ 

oo 

C e = fn(a* m a m - a*_ m a_ m ) = 0, (7) 

m=l 



C G = = Cbi + C { = 0, (8) 
2 

CBI = -^jfWCcrPcr + CgpgCsps + C a p a C S ps], (9) 

(10) 

Cbi is the classical Hamiltonian constraint of the Bianchi I model, whereas is the 
inhomogeneous term, which couples in a nontrivial manner the homogeneous sector of 
the Gowdy model with the inhomogeneities. These are encoded in the contributions 

H o = Yl H<An, H^ nt = ^ [2a* m a m + a m a_ m + a* m a*_ m ] . (11) 

These terms represent, respectively, the Hamiltonian corresponding to a free massless 
scalar field and an interaction term quadratic in the field. 

The choice of variables that we have made is the most suitable one for the subsequent 
hybrid quantization that we will perform. On the one hand, the homogeneous sector has 
been described in the Ashtekar formulation in order to prepare it for the polymeric quan- 
tization characteristic of LQC On the other hand, for the inhomogeneous sector we have 
chosen the parametrization introduced in Ref. [131 ], since it is essentially the only one 
which admits a satisfactory Fock quantization, including a unitary implementation of the 
dynamics after completing the choice of time gauge (i.e. after a complete deparametriza- 
tion). The coupling between the two sectors shows the interest (and nontriviality) of 
our approach because, owing to it, it is not straightforward that a well-defined hybrid 
quantization be viable. 



r - G 



{c a p 



3. FOCK QUANTIZATION 

We will start the construction of our quantum model with the Fock quantization of 
the inhomogeneous sector and the imposition of the constraint defined by the generator 
of ^-translations, which depends only on the inhomogeneities. In doing so, we first 
choose as complex structure the one which is naturally associated with the identification 
of {(a m , a* m )} as creation and annihilation operators, for all modes with m^O [ijj]. Using 
this complex structure, we construct the "one-particle" Hilbert space, the corresponding 
symmetric Fock space J 7 , and a quantization of the variables {(a mi a* m )} by standard 
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methods. The field dynamics is unitarily implemented in this quantization [13j. Besides, 
the quantization provides also a natural unitary implementation of the gauge group of S 1 - 
translations, since the resulting vacuum is invariant under that group. Imposing these two 
unitarity requirements, the Fock quantization turns out to be unique (up to equivalence) 

Q. 

We will denote by |{n TO }) := |..., n_ m , n m , ...) the corresponding n-particle states, 
which provide an orthonormal basis for the Fock space. Here n m < oo is the occupation 
number of the m-th mode. In these states, only a finite set of these occupation numbers 
differ from zero. The dense set spanned by them will be denoted by S. 

The generator of S^-translations, given in Eq. ([7]), can be promoted to the quantum 
operator: 

oo 

C e = ^m(aj n a m - at m a_ m ). (12) 

m=l 

The proper Fock subspace annihilated by this constraint will be called T p . A basis for it 
is provided by the n-particle states which verify the condition 



oo 



22 m ( n m - n- m ) = 0. (13) 

m=l 

Clearly, the vacuum is one of these states. 

The subspace T v is the physical Hilbert space of the deparametrized system quantized 
in Ref. jl3j ]. There, Cg was the only quantum constraint. In order to get the physical 
Hilbert space now, however, we need to impose also the quantum counterpart of the 
Hamiltonian constraint ([8]). 



4. QUANTIZATION OF THE BIANCHI I MODEL 

When the inhomogeneities vanish, the Gowdy model reduces to the vacuum Bianchi 
I model with T 3 topology. As a preliminary step before discussing the quantization of 
the complete Gowdy system, including both the homogeneous and inhomogenous sectors, 
we will focus our attention on the Bianchi I model, quantizing it polymerically. Here, 
we will summarize and revisit the analysis of Ref. [9(, where we accomplished the loop 
quantization of this particular (homogeneous) subfamily of Gowdy spacetimes. 
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4.1. Quantum representation 



In LQC the basic configuration variables are holonomies, whereas the basic momentum 
variables are fluxes. The holonomy along an edge of oriented coordinate length 27r/ij in 



the direction % is h^{c l ) = e^ liC '" ri 2lj, where Tj are the SU(2) generators proportional to 
the Pauli matrices, such that [Tj, Tj] = €ijkT k . On the other hand, the flux across a square 
S l of fiducial area (27r/i,) 2 and normal to the direction % is E(S l ) = Pi (/ij) 2 . 

The configuration algebra for each fiducial direction i is the algebra of almost periodic 
functions of c l , which is generated by the matrix elements of the holonomies M^c 1 ) = 
e fwc*_ w e will call Cylg the corresponding vector space. The kinematical Hilbert space 
Hkin is the tensor product of three copies T~t\- m (one for each fiducial direction) of the 
space L 2 (IRb, c?A*Bohr), where Rb is the Bohr compactification of the real line, and rf/iBohr 
is the normalized Haar measure on it |23|. In momentum representation, and employing 
the Dirac ket notation to denote the states A/^c 1 ), the Hilbert space 7^ in can be 
seen as the completion of the algebra Cylg with respect to the discrete inner product 

(Mi|/4> = <W<- 

On the basis states the action of the basic operators pi (associated with fluxes) 
and N^i. (associated with holonomies) reads 

Pi\/ii) = A-Kjllx/ii | A^|/Ui) = + (14) 

Here, lp\ = \J Gh is the Planck length. 

4.2. Hamiltonian constraint 

The classical Hamiltonian constraint of the Bianchi I model leads to the following 
form when expressed in terms of the basic variables of LQC: 



r - 2 

^BI — ^ 




A e A a I — == | + A e A 5 [ — = 1 + A CT A 




(15) 



where 



A, = - lim ^Msign^)^ - A/L 2/4 ). (16) 

The existence in Loop Quantum Gravity of a minimum nonzero eigenvalue for the 
area, A = 2y / 3vr7/p 1 , has been argued to imply that the limit — > is not feasible and 
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that, for each fiducial direction, there exists in fact a minimum nonzero edge length for 
the holonomy, fa. Therefore, in order to obtain the quantum counterpart of Aj in LQC, 
the above classical limit is replaced by evaluation at pi. On the other hand, to determine 
this minimum fiducial length we adopt the proposal presented in Ref. 8|, which leads to 
the condition pf\pi\ = A. Thus, the value of pi depends on that of \pi\, which is one of 
the degrees of freedom of the physical metric. Quantum mechanically, this translates into 
the operator relation 

fa 



(17) 

Taking then a suitable symmetric factor ordering in the quantum counterpart of (fl~6j) . we 
finally arrive at the operator 



A; 



4V A 



\Pi 



- A r -2w)sign(p i ) + sign^X-A/^ - A/L 



2m 



\Pi\- 



In principle, given relation (iT7|) . the operator Jv^. should produce a state-dependent 
shift fa(fa) on the basis states \fa). It is then most convenient to relabel this basis 
introducing an affine parameter such that the action of Af^ is defined to cause just 
a constant shift in the new label. Such an affine reparametrization is possible (see e.g. 
Ref. |9j), resulting in the following action of the basic operators on the states |t>j) of the 
relabeled basis: 



Pi\v 



) = 3 1/3 A sign( 



Vi) Vi 



2/3 



tffiilVi) = h + l). 



(19) 



Starting with pi, we can define the operators y\p~i\ and sign(pj) via the spectral theorem 
24|. It is then straightforward to compute the action on the above basis of the operator 
Aj, given in Eq. ([18]) . To promote the Hamiltonian constraint ([151) to an operator, we also 
need to represent the inverse of In this case, the spectral theorem is not applicable 

since zero is in the discrete spectrum of To overcome this problem, one introduces the 



quantum analog of the classical identity [2l| 




signfe 



-tr [ Ti/if (c l 



hi 



\Pi 



(20) 



2nGpi 

The result is a regularized operator whose action on the basis states turns out to be 
diagonal: 



\Vi) 



s/lh 



\vi), b(vi) 



3 7 / 



12 



,..|V3 



/87T 



k + i 



1/3 



Vi 



1 1/3 I 



(21) 
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In this way, the Hamiltonian constraint of the Bianchi I model is represented by the 
symmetric operator 



C 



BI 



T 



1 



\Ps 



A 9 A 5 



1 



A„A* 



1 



(22) 



which is densely defined on the domain Cyl s = span{|t>0, v a , v$) = \v$) ® \v a ) ® \vs)}. 



4.3. Singularity resolution 

Owing to the factor ordering chosen in our definitions, our symmetric Hamiltonian con- 
straint annihilates all the states in the basis { \v$, v a , v$)} with any of the v^s equal to zero. 
We will call them "zero volume states" , since they are eigenstates of the (homogeneous- 
volume operator V = ®i^/\pi\ with vanishing eigenvalue. Furthermore, the complement 
of the subspace spanned by these zero volume states is invariant under the action of the 
constraint operator (1221) . Therefore, the restriction of the constraint to this complement 
is well defined. We can say that, when the constraint is imposed, this complement de- 
couples from the space of zero volume states. Then, we limit our considerations to it in 
the following in order to find nontrivial solutions to the constraint. We will call Cyl 5 the 
linear span of tensor products of states |%) such that none of the v^s vanishes, whereas 
H-Km will denote the corresponding Hilbert space of nonzero volume states. 

It is worth noticing that, since zero volume states have been removed from our kine- 
matical Hilbert space, so that the kernel of all the operators pi is empty, there is no longer 
any quantum analog of the classical cosmological singularity, where some of the triads Pi 
vanish classically. In this sense, the singularity is resolved in our Bianchi I model. We 
will discuss this issue in more detail in Sect. 6, where we will also analyze the fate of the 
singularity in the Gowdy model. 



4.4. Densitized Hamiltonian constraint 



In principle, solutions to the Hamiltonian constraint do not need to be normalizable 
in the kinematical Hilbert space, which is mainly a mathematical tool introduced to con- 
struct a representation of the holonomy-flux algebra and the constraints of the model, but 
which does not take into account the dynamics. We will look for solutions in the much 
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bigger space (Cyl s )*, the algebraic dual of the dense domain of definition of the Hamil- 
tonian constraint, whose elements will be denoted by Instead of looking directly for 
solutions to the constraint (J22l) . which are difficult to determine, we find more convenient 
to densitize this constraint by means of the following bijection in the dual space (Cylg)*: 



1 

V 

where the operator that represents the inverse of the (homogeneous-) volume is 



(23) 



1" 




1 


y_ 


= ®i 


V\Pi\ 



(24) 



The transformed physical states are now annihilated by the (adjoint of the) symmetric 
densitized Hamiltonian constraint, defined as 



C 



BI 



" 1 " 


2 


" 1 " 


2 


y_ 


Cbi 


y 


7 2 



+ e e e s + e ff e a 



(25) 



where Gj is the symmetric operator 





i 


2 


1 




_V\Pi\_ 




V\Pi\ 



(26) 



Note that the operator [1 / A/bi|] -1 ^ 2 is wen defined once we have restricted the study 
to (a dense domain in) the space of nonzero volume states. It is also worth noticing 
that the above kind of densitization procedure can be similarly applied if Cyl s is re- 
placed with another dense domain for the Hamiltonian constraint such that the inverse 
(homogeneous)-volume is still a bijection in the corresponding dual. 

From Eq. (I25I) . we see that all the operators Gj are Dirac observables in the Bianchi 
I model, because they commute with the constraint. This fact simplifies enormously the 
resolution of the constraint, which can be seen as an algebraic equation in the eigenstates 
that are allowed simultaneously for the three operators Gj, one for each direction. Since 
the three operators are formally identical, the problem is reduced to a one-dimensional 
problem. With this motivation, we will now analyze the properties of Q,. 
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4-4- 1- Superselection and no-boundary description 



As one can check by direct calculation, the operator 0j is a difference operator acting 
on the basis states |t>j): 



where 



and 



%i\ Vi ) = -z-^= [f+(vi)\vi + 2) - f-(vi)\vi - 2)] 



f±(vi) = g{v,i ± 2)s±(vi)g(vi), 8±(v t ) = sign(^ ± 2) + sign(v; 



(27) 



(2f 



1 + i 



1-i 



if «i^0, 
if u = 0. 



(29) 



The presence of signs in Eq. (1281) has important consequences. Namely, it turns out that 
the function (/_(fj)) vanishes in the whole interval [—2,0] ([0,2]). As a result, in 

terms of the label Vi, the action of the difference operator Bj does not mix any of the 
following semilattices: 



C% = {He, + 2k), k = 0, 1, 2...}, e % e (0, 2]. 



(30) 



From now on, we will call TCf. the corresponding subspaces of states with support in these 
semilattices. Then, each of them provides a superselection sector. 

If we compare the classical densitized constraint [see Eqs. (jSJ) and Qj] with its quantum 
counterpart, given in Eq. ([25]). we conclude that the operator 9j is the quantum analog 
of the classical quantity dpi 21]. Hence, the Wheeler-DeWitt analog of 0j would be the 
first order differential operator Q { = i8ir'ylp l p i d pr One would expect that Q { could be 
recovered from 0j in a suitable semiclassical limit. However, its action ( 1271) indicates that 
0j is instead a second-order difference operator. This apparent conflict, nonetheless, does 
not break physical consistency. Because of the property f^(si) = 0, when one hits the 
origin in TCf., one gets a relation which constraints the data at Vi = ±£j and V{ = ±(£j+2). 
Hence, the eigenstates of 0j in Tlf. depend in fact on a single piece of initial data (as it 
would correspond to a first-order operator), given just by the projection on the Vi = ±e$ 
"slice" of the semilattice. 



14 



In this sense, besides, the Hamiltonian constraint provides a no-boundary description: 
physical states will not only have no contribution from the slice Vi = (for any i=l, 2 or 
3), which would correspond to the classical singularity, but they also get no contribution 
with the opposite sign of the label vi -namely the opposite triad orientation- without the 
need to impose any boundary condition at the initial slice. 



4-4- Spectrum and eigenf unctions 



Let us call Cyl* the linear span of the t^-states in the semilattice Cf.. Then, the 

9|. In addition, its spec- 



operator 0j, with domain Cyl E ., is essentially self-adjoint on TC Z £ 
trum has been completely characterized: it is absolutely continuous, nondegenerate, and 
coincides with the real line. Furthermore, the generalized eigenfunction with generalized 



eigenvalue Xjlpi, denoted by e 



is completely determined by the data on the initial 



slice, e A £i (±£i), as we pointed out before. One can show that the explicit expression e.g. 



of e%(vi) is 



e^ + 2M) = ]T 

O(M) 



n 



/_(e, + 2r fc + 2) 
/+(£< + 2r k + 2) 



jj -i 2 V / 3A7Z| >1 



{«!} 



Af + (e i + 2s l 



e e j(ei). (31) 



Here, O(M) denotes the set of all possible ways to move from to M by jumps of one or 
two steps. For each element in O(M), {r k } is the subset of integers followed by a jump of 
two steps, whereas {si} is the subset of integers followed by a jump of only one step. Note 
that, up to a constant phase, these complex coefficients oscillate from real to imaginary 
when Vi varies along the considered semilattice. 

After a suitable (delta-)normalization of the generalized eigenstates \e^ 6i ) (in ket nota- 
tion), the spectral resolution of the identity in the kinematical Hilbert space TCf. associated 
with 0j is given then by 



I 



± 



dA|e A £ ')(e> 



(32) 



4.5. Physical states for Bianchi I 

Let us restrict our study to any specific superselection sector determined by the three 
numbers (eg, e a , eg) and by a sign in each direction for the orientation of the triad. In 
the following, for simplicity, we will choose positive orientations. Since the Bj's are Dirac 
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observables and we know their associated resolution of the identity, it is really easy to 
find the physical Hilbert space. We can follow two strategies. On the one hand, since 
0j is essentially self-adjoint in CyLt, the constraint operator Cbi, given in Eq. ( l25|) . is 
essentially self-adjoint in the tensor product of these spaces, and we can apply the group 



averaging procedure |25l . l26j to determine the solutions to the Hamiltonian constraint 
and their Hilbert structure. On the other hand, we can also solve directly the constraint 
and determine the Hilbert structure of the solutions by choosing a complete set of real 
observables and imposing that they be represented as self-adjoint operators. 

In both cases, one concludes that the solutions to the Hamiltonian constraint have the 
form 

ip(v 6 ,v a ,v s ) = / d\ a d\ s el 9 e[X] (ve)e e £(v a )e £ £(v s )$(\ a , X s ), (33) 
Jm. 2 

with 

A,[A] = (34) 
Here, the physical states ip(X a , Xs) belong to the Hilbert space 

H B1 = L 2 (R 2 , |A CT + X s \dX a dX s ) . (35) 

4-5.1. Evolution and observables 

In the above expression for the solutions, we have eliminated the dependence on Xg, 
determined in terms of X a and A^. This is just a matter of convention: we could choose 
to fix any of these three variables in terms of the other two. In view of our choice, we can 
interpret vg as an internal time and regard physical states as evolving with respect to it 



271 ] . At a particular time Vq, the solution can be written as 

4>{v a ,v s )\ v o = dX a dX & e £ { (y (v a )e £ ^ s (v 5 )i)(X a ,X 5 )\ v o 1 (36) 



where 

^(A CT ,A,)|,o = e ^ [A] K)^(A (T ,A 5 ) (37) 
belongs to the "fg-slice" Hilbert space 

H v o = L 2 ( M 2 , |A " + Xs \ dX a dXs ) . (38) 

9 V J 
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A complete set of observables is given by the constants of motion O5 and CT , which 
act on the physical states just by multiplication by Xo-jlpi and Xsjlpi respectively, and by 
the observables at fixed time v a \v° and vs\ v %, whose action on solutions is given by 

v a \v° e i>(va,v s )\v% = v a i/j(v a ,v s )\v%, a = a,5. (39) 

On the Ug-slice Hilbert space, on the other hand, the corresponding action is 

Va\v° g ${K, X s )\ v o = J dX a (e £ {Jv a \e £ sJ^(X a , \s)\ v °, (40) 

and similarly for Vs\v°- 

The fact that the dependence of |e A V]( u 0)| on \ a and Xg varies with the value of ve 
precludes one establishing a unitary relation between the different fg-slice Hilbert spaces, 
obtained for other choices of v g instead of Vg, so that one does not get a unitary evolution 
in this emergent time. This result is a mere consequence of the nature of the polymeric 
quantization performed on the variable that plays the role of internal time 



5. HYBRID QUANTIZATION OF THE GOWDY MODEL 



Once we have represented the homogeneous sector of the Gowdy model following the 
prescriptions of LQC, and the inhomogeneous sector employing the Fock quantization, 
we are ready to construct the quantum counterpart of the Hamiltonian constraint of the 
Gowdy model, Cq, within our hybrid approach. The kinematical Hilbert space is just the 
tensor product of the kinematical Hilbert spaces for each sector, that is ?^kin <8> T . 

Taking into account Eqs. fl8])- (fTTj) . it is straightforward to represent the constraint Cq 
as an operator. In fact, in the previous section we have already represented the quantum 
analog Cpi of the homogeneous term Cpi- In addition, we have constructed the inverse 
homogeneous- volume operator, given in Eq. (}24"1) . and obtained the quantum counterpart 
Oj of the quantity dpi. Therefore, we can also promote the inhomogeneous term to 
an operator. Adopting the same kind of factor ordering used for the homogeneous term, 
we get the following symmetric operator: 



C f 



1 

V 



c F 



1 

V 



] 2 



V\Pe\ 



HL + 327T 2 



\m 



(41) 
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where the free Hamiltonian Hq and the interaction term H^ nt are normal ordered: 

#0 = Y M^An, = Y ^^ Ta ^ hm + + °L^-m)- (42) 

Similarly to what happens with the Hamiltonian constraint of Bianchi I, Cbi, the in- 
homogeneous term annihilates the subspace spanned by the zero homogeneous-volume 
states and leaves invariant its complement. Therefore, this subspace decouples from its 
complement also in the complete Gowdy model when the constraint is imposed. Thus, 
we can again restrict our study to the kinematical Hilbert space Hki n i n the homogeneous 
sector. Moreover, the Hilbert spaces ®iHf_ ®F are also superselected by the Hamiltonian 
constraint of the Gowdy model, so that for practical purposes we can further restrict 
our discussion to any of these kinematical Hilbert spaces. We continue choosing posi- 
tive orientation of the triads for simplicity and, in principle, define our operators in the 
corresponding dense set 

(g)i CylJ ® <S = span{|«0,v o .,us) ® \{n m }); v t E £+}. (43) 

As in the Bianchi I case, it is preferable to densitize the Hamiltonian constraint, in 
particular because it is then straightforward to recognize some Dirac observables, what 
facilitates the resolution of the constraint. We introduce again the bijective map ( 1231) 
[in principle in the dual (®jCyl^ <g> S)*, but the dual of any other dense set where the 
inverse homogenous- volume operator provides a bijection would be acceptable as well]. 
This leads to "transformed" physical states which are annihilated just by the densitized 
Hamiltonian constraint Cq = Cbi + Q, with Cbi given in Eq. (1221) and Q in Eq. (|4TT) . It 
is worth noting that, although the constraint operator Cq couples both the homogenous 
and inhomogeneous sectors of the Gowdy model in a highly nontrivial way, it is indeed a 
well-defined symmetric operator with domain ®jCyl^ ® S. 

In the Gowdy model, the operators and 0,5 are still Dirac observables. Nevertheless, 
the densitized Hamiltonian constraint depends now on \p$\ and l/y/\pg\, so that Qg no 
longer commutes with this constraint and fails to be a Dirac observable. It is then 
convenient to use the basis of states \v$) <S> |e^) ® |e^) for the homogeneous sector. 
In doing so, the densitized constraint becomes a difference equation in vg and can be 
regarded as an evolution equation in this parameter, what means that pg plays the role 
of internal time. As we have already commented, this was essentially the time choice 

18 



adopted to deparametrize the system in the Fock quantization of Ref. 13J. Apart from 
a global factor and some suitable rescalings, the inhomogeneous part of the constraint, 
Q, coincides in each "generalized eigenspace" of CT and 9,5 with the inhomogeneous 
Hamiltonian of the deparametrized (gauge-fixed) model [131 . The time dependence of 
that Hamiltonian becomes a dependence on pg in our case. 



5.1. Solutions to the eigenvalue equation for the Hamiltonian constraint 



Let us consider the (complex) eigenvalue equation for Cq, 



(44) 



Using the basis of states introduced above for the homogeneous sector and substituting 
the formal expansion 



vgECtf. 



dX a dX s (v e \ ® {e%\ ® (e£| ® (ipx„,\ s M\, 



(45) 



we get the solution [171, |l9] 



(:<l>x,.x,(^ + 2M)\{n m }) = (Va„,a,M Y. 



n 



/_(e e + 2r fc + 2) 
f+(e e + 2r k + 2) 



O(M) 

x v[Y[H^ee + 2 Sl ,X a ,X s )\ \{n m }), 



(46) 



where O(M), {r^} and {s;} have the same meaning as in Sec. (14 4.4 4.4.21) . Besides, the 
symbol V denotes path ordering, and Hj-(vg, X a , Xs) has the form: 



H^(v e , X ff , X s 



27r(X a + X s )f + (v e ) 



P + 2X a X s 



(X a + X 



S) 1,2 
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(47) 



5.2. Physical states and observables 



The operator Cq is essentially self-adjoint if and only if there is no normalizable solution 
of the form (T46l) when p = ±z. In principle, one can make use in that case of the 



group averaging method 



25 



26] to construct the physical Hilbert space. In practice, 
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nevertheless, to apply this method we should first determine and adopt an invariant 
domain for the constraint operator Cq, a task which seems difficult to accomplish given 
the complexity of the action of the constraint and the infinite number of degrees of freedom 
of the system. Therefore, instead of following that approach, we will directly employ our 
knowledge of the formal solutions to the densitized Hamiltonian constraint, deduced in 
the previous subsection, and determine a physical Hilbert structure for them by requiring 
the self-adjointness of a complete set of real observables. 

Solutions to the constraint are formally given by expression (|46p with p = 0. As we 
can see, the initial data (^>x a ,x s ( £ e)\ completely determines the formal solution, so that 
we can identify the latter with the data on the initial section vq = eg. To construct the 
physical Hilbert space, we then simply provide the vector space of initial data with a 
Hilbert structure. The desired inner product can be fixed by choosing a complete set 
of real classical observables and requiring that their quantum analogs, which act on the 
initial data, be self-adjoint operators. Such a complete set is provided by the observables 
introduced for Bianchi I in Sec. (144.5 4.5TTI) (with a trivial action on the inhomogeneous 
sector) and by a complete set acting on the inhomogeneous modes, which can be, for 
instance, the set of operators that represent the Fourier sine and cosine coefficients of the 
nonzero modes. Up to irrelevant constant real factors, these operators are 

{(aw + ajj ± (a_ m + al m ), i[(a m - ajj ± (a_ m - at m )]; m G N + | . (48) 

Actually, they are self-adjoint operators in the standard Fock space T . Therefore, we 
conclude that the Hilbert space of initial data picked up by our conditions is, up to 
equivalence, L 2 (IR 2 , (A^ + \$\d\ a d\s) £g> T. 

Finally, to obtain the true physical Hilbert space of the model, we still have to impose 
the ^-symmetry generated by the constraint Cg, which commutes with the Hamiltonian 
constraint. This symmetry is encoded in the condition (|T3l) . Taking this into account, we 
arrive in the end to the physical Hilbert space 

H G = L 2 (R 2 , | A CT + Xs\dX a dX s ) ® T v . (49) 
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6. CONCLUDING REMARKS 



We have carried out a thorough quantization of the Gowdy spacetimes with the spatial 
section of a three-torus and linear polarization by combining the loop quantization of the 
homogenous Bianchi I cosmology (with compact sections) [9( with the Fock quantization 
of the inhomogeneities in a totally deparametrized Gowdy system 13|. In this way, we 
have constructed a hybrid quantum model for this family of cosmological spacetimes in 
vacuo which incorporates the presumably most relevant effects of the quantization of 
the geometry (at least in scenarios that can be considered close to homogeneity), while 
allowing the treatment of an infinite number of degrees of freedom. 

Even though the Hamiltonian constraint couples in a nontrivial manner the polymeric 
quantum homogeneous sector with the inhomogeneous sector, quantized with standard 
Fock methods, we have been able to find the formal solutions to this constraint. The 
constraint can be regarded as an evolution equation in an emergent time, and the explicit 
expression of the solutions shows that they are in fact completely determined by their 
data at an initial value of this time. 

The emergent time corresponds to one of the variables that have been quantized with 
loop techniques. Because of this fact, a naive and straightforward definition of the evo- 
lution with respect to it is not implemented as a unitary transformation in the quantum 
theory. We note, in this sense, that the results discussed at the end of Sec. (14 4.5 4.5.11) 
for Bianchi I may be extrapolated (at least in principle) to the case of the hybrid Gowdy 
model, because this model admits observables which are just the tensor product of Bianchi 
I observables and the identity on the inhomogeneous sector. 

Owing to the polymeric quantization carried out in the homogeneous sector, we have 
been able to decouple the space of states in the kernel of any of the triad operators and 
rigorously remove it from our quantization. Therefore, there is no longer a quantum analog 
of the (generic) classical cosmological singularity in the system. This result, originally 
obtained for the Bianchi I model, extends also to the Gowdy model. Hence, we conclude 
that the loop quantization of the considered zero modes suffice to avoid the singularity in 
this inhomogeneous model. This robust result is different and simpler than other suggested 
possibilities for avoiding cosmological singularities in inhomogeneous scenarios, like e.g. by 
appealing to the BKL (Belinsky, Khalatnikov, and Lifshitz) conjecture 28j. In particular, 
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the persistence of the resolution mechanism found in homogenous LQC is a global result, in 
the sense that one does not need to analyze the approach to the singularity independently 
at each point of the corresponding spatial section. Furthermore, no particular boundary 
condition has to be imposed in the construction of the quantum states in order to avoid 
the singularity and prevent the emergence of contributions with different triad orientation. 
From this perspective, one can say that the developed formalism provides a no boundary 
description. 

In addition to all the above, we have been able to complete the hybrid quantization by 
determining the Hilbert space of physical states and providing a complete set of observ- 
ables acting on it. Remarkably, the physical Hilbert space which results of this hybrid 
quantization is (equivalent to) the tensor product of the physical Hilbert space for Bianchi 
I in LQC and the Fock space which describes the quantization of the inhomogeneities in 
the totally deparametrized system [Lj]. Therefore, we can say that one recovers the stan- 
dard quantum field theory for the inhomogeneities on a polymerically quantized Bianchi 
I background. Finally, let us mention that there are other interesting issues for future 
investigation for which our approach may be specially appropriate, like e.g. the analysis 
of the semiclassical behavior of physical states, or the implementation of perturbative 
approaches which deal with the inhomogeneities as perturbations around the Bianchi I 
background. 
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